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A covariant quantization scheme employing reducible representations of canonical commutation 
relations with positive-definite metric and Hermitian four-potentials (an alternative to the Gupta- 
Bleuler method) is tested on the example of quantum electromagnetic fields produced by a classical 
current. The Heisenberg dynamics can be consistently formulated since the fields are given by 
operators and not operator-valued distributions. The scheme involves a Hamiltonian whose free part 
is modified but the minimal-coupling interaction is the standard one. Solving Heisenberg equations 
of motion under the assumption that the fields are free for times to = ±oo we arrive at retarded and 
advanced solutions. Once we have these solutions we can deduce the form of evolution of retarded 
and advanced fields between two arbitrary finite times. The appropriate unitary evolution operators 
are found and their generators are computed. Now the generators involve the same free part as 
before, but the interaction term turns out to be modified. For a pointlike charge localized on a world- 
line z a (t) we find the interaction term of the form — qA(z(t)) ■ v(t) — q J dz- E where the integration 
is along those parts of the charge world-line where the charge velocity is nonzero. There is no self- 
energy contribution. Next we compute photon statistics. Poisson statistics naturally results and 
infrared divergence can be avoided even for pointlike sources. Classical fields produced by classical 
sources can be obtained if one computes coherent-state averages of Heisenberg-picture operators. 
It is shown that the new form of representation automatically smears out pointlike currents. We 
discuss in detail Poincare covariance of the theory and the role of Bogoliubov transformations for the 
issue of gauge invariance. The representation we employ is parametrized by a number that is related 
to Renyi's a. It is shown that the "Shannon limit" a — ► 1 plays here a role of a correspondence 
principle with the standard regularized formalism. 

PACS numbers: 03.70.+k, 41.20.Jb, 42.50.-p 

I. INTRODUCTION 

Dynamics of observables is given in quantum mechanics by Heisenberg equations of motion. In quantum field 
theory one does not really work with Heisenberg equations but prefers the S-matrix formalism. The situation is 
caused mainly by the fact that quantum fields are operator valued distributions and their products taken at the same 
point in space-time are meaningless. Formal solutions, when inserted back into Heisenberg equations, typically lead 
to mathematical absurdities due to divergences. The difficulties are less visible and thus easier to live with at the 
S'-matrix level. Another reason for the popularity of the S'-matrix formalism is the question of gauge invariance. There 
are various 'proofs' in the literature that the S-matrix formalism of quantum electrodynamics (QED) is manifestly 
gauge invariant. However, as discussed in detail in 1], all the 'easy one-line proofs', including the celebrated Feynman 
one |2|, are false if divergent expressions are encountered. 

As often stressed by Dirac, the removal-of-divergences techniques practically mean a departure from the Heisenberg 
equation. Dirac regarded this as a serious drawback of quantum field theory. In his last published lecture entitled 
'The requirements of fundamental physical theory', and given to a gathering of Nobel Laureates at Lindau on July 1, 
1982, he said: 

"I feel that we have to insist on the validity of this Heisenberg equation. This is the whole basis of quantum theory. 
We have got to hold onto it whatever we do, and if the equation gives results which are not correct it means that we 
are using the wrong Hamiltonian. This is the point I want to emphasise (...). Heisenberg originally formulated these 
equations with the dynamical variables appearing as matrices. You can generalise this very much by allowing more 
general kinds of quantities for your dynamical variables. They can be any algebraic quantities such that you do not 
in general have commutative multiplication (...). Some day people will find the correct Hamiltonian and there will be 
some new degrees of freedom, something which we cannot understand according to classical ideas, playing a role in 
the foundations of quantum mechanics." 

The talk appeared in a published form in 1984, the year of Dirac's death. Also in 1984 Dirac published his last 
paper [4|, entitled 'The future of atomic physics' where he criticised, again from the perspective of the divergences, 
the dogma of irreducibility of representations of elementary quantum symmetries. Taken together, the two papers 
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form a kind of scientific testament of this great physicist. 

Quite recently, in a series of papers [EIEQ) one of us was advocating the idea that the occurrence of the divergences 
may be related to the fact that the standard quantization scheme is based on irreducible representations of canonincal 
(anti-)commutation relations. It was argued that there are physical reasons for the use of certain reducible represen- 
tations of CCR (in jalfl) an d CAR (in pfl). There is no problem with multiplying fields at the same space-time point 
because the new representation leads to fields that are operators and not operator- valued distributions. Preliminary 
analysis of interactions with charges discussed in these papers (two-level atoms p|, classical current [(J, reducibly 
quantized fermionic fields 0) was always pointing into the conclusion that the reducible representation may indeed 
remove divergences, although definite statements would have been far premature at that stage. 

Particularly intriguinging seems the possibility, discussed in 0, Q , that the experiments on Rabi oscillations per- 
formed by the Paris group |10| may be more consistently interpreted in terms of fields quantized by means of the 
reducible representations. If this were the case, it would mean that the fundamental question which representation 
is more physical would be within the reach of present-day optical experiments. But then we have to go further with 
the analysis of the new representation. In particular, we have to control the relativistic and gauge properties of the 
formalism, a fact that justifies the origin of the present paper. 

The goal of the present work is to apply a covariant analogue of the representation introduced in [6| to the problem 
of Heisenberg-picture evolution of electromagnetic fields interacting with a classical current. The issue of infinities is 
not our main concern here since any reliable discussion of ultraviolet divergences would require quantized currents, 
but various automatic regularizations due to our choice of representation can also be observed. In the present paper 
we are more interested in the problem of Poincare covariance vs. gauge invariance, positive-definiteness of the scalar 
product, and unitarity of evolution. Interaction with classical currents does not grasp all the possible subtleties of 
QED, but it allows for exact solutions and thus is a natural playground for testing any new quantization paradigm. 

The results are quite promising. As opposed to standard quantization schemes it makes in our representation perfect 
sense to speak of the Heisenberg-picture evolution of field operators. We begin with solving the equations under the 
assumption that the fields are free at timelike infinities. We make this condition precise by first taking solutions that 
are free at an arbitrary to, and afterwards taking the limit to — > ±oo. The results of this limiting procedure are, 
respectively, retarded and advanced solutions of the Maxwell equations. Having the solutions for all times, we can find 
a unitary map W±(t,t\) that maps retarded/advanced fields at time t\ into retarded/advanced fields at time t. Once 
we have this unitary map, we can compute its generator. We show that it differs from the usual minimal-coupling 
Hamiltonian. 

We study in more detail the solutions corresponding to pointlike classical sources. Taking a relatively general form 
of the charge world-line we compute the explicit form of the modified interaction Hamiltonian. The interaction term 
does not contain any self-energy contribution but we find in addition a term that describes the work performed by 
the charge against free electric field. We then show what kind of a classical field one arrives at if one computes 
coherent-state averages of our field operators. The resulting field looks as if it were produced by an extended current, 
even though the current is pointlike. We construct the evolution operator and consider its associated S'-matrix. No 
infrared divergence occurs and there is no problem with computing photon statistics even for pointlike sources. 

What is interesting, the different form of representation of CCR we work with implies also a slightly modified 
form of the free part of the Hamiltonian. The modification is subtle and is related to some additional degree of 
freedom characterized by a parameter TV. The role of N for Rabi oscillations was discussed in detail in [8( with 
the conclusion that finite N are compatible with experiment. The "thermodynamic limit" N — > oo plays a role of 
a correspondence principle that maps the new theory into a regularized form of the standard one. Effectively, the 
limiting results are equivalent to the standard ones but with a cut-off. It has to be stressed, however, that there is 
no cut-off in the Hamiltonian, and the regularization is a result of the special form of the vacuum one employs in the 
new representation. We also show that the limit N — > oo is equivalent to Renyi's limit a — » 1. 

The modification of the formalism we discuss seems to satisfy standards which are not that far from the requirements 
proposed by Dirac: We have new types of noncommuting dynamical variables, they involve certain new degrees 
of freedom that have no counterpart in classical electrodynamics, and their dynamics is governed by a modified 
Hamiltonian. The dynamics is given by Heisenberg equations but their solutions simultaneously satisfy Maxwell's 
equations with a regularized current. The regularization is not introduced ad hoc, but follows from the quantum 
structure. As opposed to standard quantization schemes the reducibly quantized fields are less singular than their 
classical counterparts. One can also think of this effect as an example of David Finkelstein's idea of 'quantization as 
regularization' \V\\. 

The paper is organized as follows. We choose the Penrose-Rindler spinor formalism since it naturally leads to 
null tetrads that are implicitly present in Lorenz-gauge potentials, and play a role of polarization vectors. We begin 
with clarifying links between the nonuniqucness of Lorenz-gauge classical 4-potentials and equivalence classes of spin 
frames associated with null 4-momenta. The spin frames are later used in construction of two types of momentum- 
dependent tetrads (associated with circular and linear polarizations of spin-1 fields and 'longitudinal' and 'timelike' 
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polarizations of two additional scalar fields). As an intermediate step towards quantization we explain in Section V 
how a change of spin frame within an equivalence class is related to a Bogoliubov transformation of creation and 
annihilation operators. For classical fields we do not have to worry about changes within equivalence classes since 
anyway the result is a gauge transformation. However, when we quantize the 4-potential the changes of spin frames 
become visible in transformations of the two 'timelike' and 'longitudinal' fields and the formalism becomes ambiguous 
unless we compensate this modification by a unitary transformation of annihilation and creation operators — this is 
where we need a Bogoliubov transformation. In Section VI we introduce the reducible representation, parametrized by 
a natural number TV, of canonical commutation relations and, in Section VII, we briefly explain why in thermodynamic 
limit N — > oo the representation will automatically introduce cut-offs at the level of averages. In Section VII we show 
how to quantize the 4-potcntial in a manifestly covariant way and without indefinite metric or non-Hermitian field 
operators. In Section IX we discuss the analogue of the Jordan-Pauli function and show that distributions typical 
of the standard formalism are here replaced by non-singular objects. Section X discusses Poincare covariance of the 
theory Sections XI and XII investigate various subtleties of the Lorenz condition. In Section XIII we introduce 
vacuum, multi-photon, and coherent states. The problem of fields produced by classical currents is discussed in 
Section XIV. We solve Heisenberg equations under the assumption that at to the fields are free, and we conclude that 
the assumption makes physical sense only if to is moved to ±00. In these limits the solutions of Heisenberg equations 
are equivalent to retarded and advanced fields, but the effective current that occurs in the source term is not a 'c- 
number' but an operator that nevertheless commutes with all field operators (the operator can be nontrivial because 
the representation is reducible). Having these solutions we are in position to find an explicit unitary transformation 
that shifts a retarded/advanced solution in time by a finite At. We find both the transformation and its generator 
in Section XVI, and in Section XVII we check the result on a pointlike charge. In this way we have systematically 
derived the Hamiltonian that allows to solve the problem with retarded or advanced initial condition at t = instead 
of that with free field at ±00. In Section XVIII we address the issue of S'-matrix and the associated photon statistics. 
Finally, having the operator solutions we compute their coherent-state averages and discuss the implications of the 
formalism for classical electrodynamics. 

II. EQUIVALENCE CLASSES OF SPIN-FRAMES AND CLASSICAL GAUGE FREEDOM 

The 4-momentum k a = k a (k) = (|fe|,fe) of a massless particle can be written in spinor notation as k a — 
TT A (k)Tt A (k), where n A (k) is a spinor field defined by k a up to a phase factor. For any w (k) there exists another 
spinor oj A (k) satisfying the spin- frame condition uj A (k)ir A (k) = 1. Given ir A (k) we cannot find a unique u A (k), since 
for any function 4>{k) the new field 

Q A (k) =uj A {k)+4>{k)n A (k) (1) 

also satisfies Cj A (k)ir A (k) = 1. This leads to the equivalence relation: uj A (k) ~ uj A (k) iff ui A (k) — ui A (k) is proportional 
to Tr A (k). 

Free classical electromagnetic fields are related to ir A (k) by 

F ab (x) = d a A b (x) - d b A a (x) = -J dT(k)n A (k)ir B (k)e A , B > (a(k, -)e"* a: + a(k, +)e lk ' x ) + c.c. (2) 

where dT(k) is the invariant measure on the light-cone. The 4-potential in a Lorenz |l3j | gauge can be taken in the 
form (cf. [lj ) 

A a (x) = iJdT(k)e-' lk - x (u JA (k)Tf A ,(k)a(k,+)+7T A {k)u, AI (k)a(k,-)) + c.c. (3) 
Now, if we replace tu A (k) by uj A (k) belonging to the same equivalence class, i.e. satisfying Q), then 

A a (x) » A a (x) = A a (x)-d a $(x) (4) 

where 

$(x) = JdT(k)<j)(k)(a(k,+)e- lk - x + a(k,-)e lk - x ) + c.c. (5) 

is a solution of OQ(x) = 0. It follows that the equivalence class of spin-frames corresponds to an equivalence class of 
Lorenz-gauge potentials. 
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III. TETRADS ASSOCIATED WITH SPIN-FRAMES 

Following 01 we introduce two types of tetrads associated with the spin- frames. The null tetrad employs the 
4-momentum as one of its elements: 

k a (k) = n A {k)Tr A ,{k) (6) 

uj a (k) = uj A (k)uj A ,(k) (7) 

m a (k) = u) A (k)n A ,(k) (8) 

fa a (k) = n A (k)Q A ,(k). (9) 

Complex 4- vectors m a (k), fh a (k) occur in Q and play the role of circular polarization vectors. The Lorenz condition 
satisfied by J3J follows from the transversality property k a (k)m a (k) — k a (k)fh a (k) = 0. The null tetrad is related to 
the Minkowski-space metric tensor of signature (+, — , — , — ) in the standard way [l2| 

g ab = k a (k)ui b (k) + uj a {k)k b {k) - m a {k)m b (k) - m b (k)fh a (k). (10) 

This formula is independent of the choice of the representative uj A (k) of a given equivalence class. 
The tetrad defined by 

x a (k) = -^(m a (k) + m a (k)) (11) 

Va{k) = ^K(fe)-m a (fe)) (12) 

z a (k) = ^(oj a (k) - k a (k)) (13) 

t a (k) = -t={uj a (k) + k a {k)) (14) 

is, in the terminology of a restricted Minkowski tetrad, satisfying 

9ab = t a {k)t b (k) - x a (k)x b (k) - y a (k)y b (k) - z a (k)z b (k) (15) 
The potential now can be written as 

A a (x) = i J , dT(k)(x a (k)a{k,l) + y a {k)a{k 7 2)) e~ ik x + c.c. (16) 
The link between the two types of amplitudes 

a(fc,±) = ^(a(fc,l)±itt(fc,2)) 
is analogous to this between circular and linear polarizations. The link is not accidental. 



IV. TRANSFORMATION PROPERTIES OF SPIN-FRAMES 



The transformation properties of spin-frames we shall discuss below do not depend on their explicit realization. Let 
us denote by Afc the spacelike part of the 4- vector A a b k b (k), and by A A B the unprimed SL(2,C) matrix corresponding 
to A a b € SO (1,3). The spinor-field transformation 

7T A (fc) ^ Air A {k)=A A B ir B {A- 1 k) (18) 

implies that k a (k) — Air A (k)Air A /(k) and, hence, 

Air A (k) = e- m(A ^7T A (k) (19) 

The phase factor 

e -i0(A,fc) = A AB w A (fe)7r B (A- 1 fe) (20) 

is the one occuring in the unitary spin-1/2 zero-mass representation of the (covering space of the) Poincare group, 
and does not depend on the choice of the representative uj A (k). The angle 6 (A, k) is known as the Wigner phase. 
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An analogously defined 

Aoj A (k) = A A B uj B {A- 1 k) (21) 

satisfies 

Auj A (k)Air A (k) = 1 (22) 

and thus 

Au A (k) = e te ^(iu A (k)+^k)TT A (k)) (23) 
= e l6(A ^u A (k) (24) 

with some 0(fc). The new oi) A (k) belongs to the equivalence class of u A (k). The gauge transformation u A (k) i— » u)^(fe) 
is in general nontrivial and depends on the explicit form of the spin-frame one works with. This is the reason why 
4-potentials are not 4- vector fields: A Lorentz transformation changes Lorenz gauges within the equivalence class of 
spin-frames. In order to guarantee independence of physical quantities of particular representatives one employs the 
result described in the next Section. 



V. EQUIVALENCE CLASSES OF SPIN FRAMES AND BOGOLIUBOV TRANSFORMATION 



Of crucial importance for the formalism we shall develop below is the link between the transformation and a 
transformation of a Bogoliubov type, i.e. the one mixing creation and annihilation operators in a way that preserves 
canonical commutation relations. 

To begin with, consider four annihilation operators aj satisfying [aj,aj,] = Sjfl and the operator 



V a = x a a x + y a a 2 + z a a 3 + t a a\. 
The change of spin-frame by u> A i— ► lo a + l^je -19 ^ is represented at the tetrad level by 



(25) 



/ ta \ 




( t a 








Va 




Va 


\ Z a J 


\ Z a 



/ 1 + H72 |0| cose |0| sine -|0| 2 /2 

|0| cos 9 1 - 10| cos 

|0| sine 1 -|0| sing 

\ |0|72 |0| cose |0| sine 1 - |0|72 




(26) 



The matrix in l|26|l is a special Lorentz transformation in Minkowski space with metric tensor diag(l, — 1, — 1, — 1) 
Let us now introduce four new operators bj defined implicitly by 



V a 



t a al + x a a x + y a a 2 + z a a 3 



x a bx + y a b 2 + z a b 3 



(27) 



The explicit form 



b\ = a\ — |0| cosea3 — |0| coseap 

62 = a,2 — |0| sinea3 — |0| sineap 

b 3 = |0|coseai -I- |0|sinea 2 + (1 - i|0| 2 )a 3 - i|0| 2 aj, 

b = -|0|cosea t 1 -|0|sine4 + i|0| 2 4 + (l + i|0| 2 )a o 

implies that = Sjfl and, hence, there exists a unitary Bogoliubov- type transformation B satisfying 

bj = B^cijB. In order to explicitly construct B 0] we introduce the representation of the Lie algebra so(l, 3), 

(28) 
(29) 
(30) 
(31) 
(32) 
(33) 



Jl 


= i{a\a 2 


- ala 3 ), 


■h 


= i{a\a 3 


- 4 a i)> 


J 3 


= i{a\ai 


- a\a 2 ), 


Ki 


= *(«o a i 


- a ai), 


K 2 


= i( a o a l 


- a a 2 ), 


K 3 


= *( a o4 


- aoa 3 ), 
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and their combinations that form a representation of e(2), 



Li 


= K x + J 2 


= i(, 


L 2 


= Ki-Jx 


= i(, 


L 3 


= h, 




[Li,L 3 ] 


= -iL 2 , 




[L 2 ,L 3 ] 


= iLi, 




[Li,L 2 ] 


= 0. 





■ 1 "T - a ai + a\a 3 — a' 3 ax), (34) 

'2 



■'■■>- a a 2 - a\a 2 + a\a 3 ), (35) 



Then 



B = B(<j>) = e i\4>\(Llcos8+L 2 sin8)^ ^ 



(36) 
(37) 
(38) 
(39) 



as can be verified by a straightforward computation. Obviously, [B^x), B(<fi 2 )] = for any 1; cf> 2 , since the map 
4> i ^ B(0) is a representation of translations. 



VI. CONSTRUCTION OF THE REDUCIBLE REPRESENTATION OF CCR 

The construction is analogous to the one introduced in [(|. The modification with respect to is that here we 
introduce four types of annihilation operators, and not just two corresponding to the polarization degrees of freedom. 

One begins with four operators, ao, a 2l a 2 , a 3 , satisfying commutation relations typical of an irreducible representa- 
tion of CCR: [oj,o!-/] = Sjj'l. Let |0) denote their common vacuum, i.e. <Xj\0) — 0. Now take the kets \k) normalized 
with respect to the light-cone delta function 

(k\k') = S r (k,k') = (27r) 3 2|fc|(5 (3) (fc-fc'). (41) 

What we call the N = 1 (or 1-oscillator) representation of CCR acts in the Hilbert space Ti spanned by kets of the 
form 

| t .no,n„n 2 ,n3) = W« (4) ~'° i n 4 r! 4 '" , |°)- 

V no ] -ni\n 2 \n 3 \ 

Physically, TL may be regarded as representing the space of states of a single four-dimensional oscillator. The 1- 
oscillator representation is defined by 

a{k,j) = \k){k\®a j . (42) 

This representation is reducible since the commutator 

[a(k,j),a(k',j^} = 6 jf 5 T (k,k')\k){k\®l (43) 

involves at the right-hand-side the operator-valued distribution J(fe) = |fe)(fe| ® 1 belonging to the center of the 
algebra, [a(k, j), I(k')] = [I(k), a(k' , j'y] = 0, for all fc, k' , j, j'. The I(k) form a resolution of unity 

JdT(k)I{k) = JdT(k)\k)(k\(g>l = I. (44) 

Here I is the identity operator in TL. 

For arbitrary N the representation is constructed as follows. Let TL be the representation space of the N = 1 
representation. Define 

TL = TL® -y®H (45) 

N 

and let A be an arbitrary operator defined for N — 1. Let 

A {n) = I® ■ ■ ■ ® I® A® I® ■ ■ ■ ® I. (46) 

V v ' V v ' 

n— 1 N—n 
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The N oscillator extension of a(k,j) is defined by 



a(k,j) = -kEn=Mk,j) {n) (47) 



and satisfies the reducible representation 

\a(k,j),a(k',j^} = 5 jf S r (k,k')I(k) (48) 

where 

L(k) = £E1i'(*) (t °- (49) 

As before we find the resolution of unity 

JdT{k)I(k) = I (50) 

where J is the identity operator in H. 

VII. THERMODYNAMIC LIMIT AND QUANTUM LAW OF LARGE NUMBERS 

The asymptotic properties of the formalism for N — > oo can be anticipated already at this stage if one reco gniz es in 
the formula 149f) the frequency operator employed in the analysis of quantum laws of large numbers |25l l26l l27l l28| . 
Indeed, let Pg be a projector. The frequency operator corresponding to the random variable (proposition) represented 
by Pg reads 

Pe,N = j?ZtiPt\ (51) 

Let Pe\0) = \6) and = J2 e ipg\8) be a state. Let = <8 • • • <8 \ip) (N times). Then the following weak law of 
large numbers holds true 

hm \\(Pe, N n)-\^\ 2m W\\=0 (52) 

for m = 1,2,3,.... The weak law states that effectively, for N — > oo, the frequency operators act by multiplication 
of ./V-copy states by appropriate probabilities. This is essentially why also in our field theory we will find averages 
where, in the limit ./V — > oo, the operators £(fe) will be replaced by their corresponding probabilities Z{k) associated 
with the choice of the vacuum state. 

VIII. QUANTIZATION OF THE POTENTIAL 

The potential operator at the level of the iV-oscillator representation reads 

A a (x) = iJdT(k)(x a (k)a(k,l)+y a (k)a(k,2) + z a (k)a(k,3)+t a (k)a(k,0y)e- lk - x + h.c. (53) 
= i J dT(k)(m a (k)a(k, +) + jfi a {k)a(k, -) + z a (k)a(k, 3) + t a (k)a(k, 0) t )e- lfe ' a; + h.c. (54) 

It is better to think of (|53|l and l|54|) as operators representing a system quantized at the N = 1 level, and then 
extended to arbitrary N by 

A a {x) - jzYLiM*)™- (55) 

From such a perspective it is easier to understand the structure of generators of the Poincare group and other 
observables. 

Let us note that the potential is Hermitian and the Hilbert space 7£ involves a positive-definite scalar product. A 
change Q of spin-frame can be compensated by an iV-oscillator Bogoliubov transformation P>_ = B ® • • • eg) B, where 
B is of the type discussed in Section Ivl 

The commutator of fields taken at arbitrary space-time points 

\Aa(x),My)} = *9abD{x-y) (56) 
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involves the operator analogue of the Jordan-Pauli function 

D(x) = iJdT(k)L{k)(e~ lk x - e* x ). (57) 

The correct signature of the metric tensor in 1(56(1 comes from the Bogoliubov-type structure of the positive-frequency 
part of il a (x), i.e. the combination of annihilation and creation operators. If one had replaced a by ao one would 
have been forced to depart either from positi vity of the scalar product or unitarity of evolution. A reducible version of 
such a (Gupta-Bleuler) formalism is possible |l7j . but one can show that contradictions with probability interpretation 
of the theory would necessary occur (for a brief discussion of the problem cf. Section IXIXfl . 

IX. JORDAN-PAULI OPERATOR: THE ROOTS OF REGULARIZATION 

To understand why the formalism we construct is less singular than the one based on irreducible representations 
it is instructive to take a closer look at 1(57(1 . In the first place, formula 1(57(1 is typical of all the representations of 
CCR, reducible or irreducible, the differences boiling down to different explicit forms of the central element Uk)- The 
standard Pauli- Jordan function corresponds to representations where /(fc) equals the identity. In our representation 
we can write (cf. Eq. (|^|) ~) 

D(x) = D i+) (x) +D ( -\x) (58) 
D {±) (x) - ±ifdT(k)I(k)e^ ik - x = ^E^=i-D (±) (^) (n) - (59) 
The operator whose TV-oscillator extensions occur in (|59|) reads explicitly 

D&>(x) = ±iJdT(k)\k)(k\e^ k - x (E)l = ±ie Tfk - x (E)lk a = JdT{k)k a \k)(k\. (60) 

As we can see, the operators D^\x) are unitary representations of 4-translations, and their generators are given by 
k a . In particular, 

£> (±) (0) = ±il, D {±) {0) = ±il (61) 

Quantization in terms of our reducible representation replaces distributions by unitary operators. This is the main 
difference with respect to the schemes based on regularizations of distributions 0, 0, |2(], El HE HE H-l • I n our 
approach there are no cut-off functions in Heisenberg-picture operators. They appear effectively at the level of averages 
and are due to the properties of states (e.g. compare the operator 1)135(1 . involving the frequency operator I_{k) and 
no cut-off, with the average l|171f) . involving the probability Z(k) and the cut-off function x(^))- As one of the 
consequences, spectra of Hamiltonians occurring in reducibly quantized theories will not depend on cut-offs. 

These facts show that the regularization occuring in our approach is of an entirely different origin than the ones we 
know from quantizations based on irreducible representations. 

X. POINCARE COVARIANCE OF FREE FIELDS 

The Poincare transformations will be taken in the form 

a(k,±) f > e ±2ie ( A <Ve ik -Va(A- 1 k,±) = u{ jy a(k,±)U A , y (62) 
a(fe,3) i ^ e ife ' !/ a(A- :L fe,3) = [/] Vi!/ a(fe,3)f/ Ai!/ (63) 
a(fc,0) t i ^ e lfc - y a(A- 1 fe,0) t = U\. y a{k^U Ay (64) 

where 9(A, k) is the Wigner phase. Transformation (|62|l is the unitary spin-1 massless representation of the Poincare 
group. Transformations l(63|l . 1(64(1 imply that the additional two fields are spin-0 and massless. Similarly to @ we 
reduce the construction to the problem of finding U\, y satisfying 

U-A, y = Ua, v ® ■ ■ • ® Ua, v/ (65) 

V 

N 
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A. Four-translations 

The 4-momentum for N = 1 reads 

P a = JdT(k)k a \k)(k\ <g) (a\ai + ala 2 + ala 3 - ala ) 



JdT(k)k a \k)(k\ (a\ ai +ala 2 ) + JdT(k)k a J(k) : (66) 



P i p ii 



with 



One immediately verifies that 



implying 



J(k) = \k) (k\ <g> (4a 3 - 4a ) (67) 



e iP ' x a(fc, ±) e -* F ' x = a(fe, ±)e _te - fc (68) 
e^a^, 3)e~* F ' x = a(fc, 3)e~ ra ' fc (69) 
e lP - x a{k^e- lP - x = a(k, OjV*-* (70) 



lA,yAa{x)U-x,y = Mx-y). (71) 



The operator J(k) will later reappear as the generator of rotations in the group E(2) associated with the 4-potential. 
The part P£ is identical to the 4-momentum operator introduced in [fj. The 4-momentum for arbitrary N reads 

P a = ZLiPP (72) 
= Ei + Jdr(k)k a J(k) (73) 

Only for N — 1 the expression coincides with the generator found by standard Noether formulas (cf. the discussion 
of this point in [5|). The form (|72|l is characteristic of a 4-momentum of N non-interacting particles. These particles 
(four-dimensional oscillators) have no counterpart in classical electrodynamics. It should be stressed that these are 
not the oscillators of Heisenberg, Born, and Jordan |29j | since there is no relationship between N, which is finite, and 
the number of different frequencies, which is infinite. 

B. Rotations and boosts 

To find an analogous representation of 

a(k,±) i — ► e ±2ie ^a(A- 1 k,±) = u} Afi a(k,±)U Ai0 (74) 
o(fc,3) h-> a(A _1 fe, 3) = ZZl i0 a(fc, 3)C7 Ai0 (75) 
a(fc,0) f h-> a(A- 1 fe,0) t = U{, o a(k,0)^U KO (76) 

we take the same definition as in Q, i.e. 

U A ,o = cxp(E s=± 2 lS /dr(fc)0(A,fc)|fc)(fc|®at as )(/ ( ir(p)|p)(A- 1 p|®l). (77) 
Taking into account the properties of spin-frames and tetrads one verifies that 

UloAtWAfi = ^JdT(k)e~^ lx 

x (m Q (Afc)e 2ie(A ' Afc) a(fc, +) + m (Afe)e- 2ie ( A > Afe) a(fc, -) + 5 a (Afc)a(fc, 3) + t a (Ak)a(k, 0) f ) + h.c. 
= A^A" 1 *) = A a fc S(A)tA h (A- 1 a ;)B(A) (78) 
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where B_( A) compensates the change of gauge caused by (|24ll . One can costruct B_(A) by first finding an appropriate 
-B(A) of the form analog ous to gnjl, and then defining 5(A) = B(A)< g > N . It is more elegant to assume that U_ A 
is accompanied by a redefinition of vacuum (see below) \Q) i— > S(A)^|0). Then the transformation of the potential 
becomes effectively 

ui, y Aa{x)U K , y = Aa'A^A-^x - y)) (79) 

i.e. that of a 4-vector field. 

A still simpler way is to assume the transformation rule 

A A B uj B (A- 1 k) = e 40(A < fc W(fc) (80) 

i.e. to replace (I24|l by the form typical of the entire equivalence class. Then one can put B_(A) = L 

The operators ocuring at right-hand-sides of field commutators transform as translation invariant scalar fields 

U{ ty L(k)U Ajy = /(A- 1 *) (81) 
Ui y D(x)U Aty = D{A- l x). (82) 

The representation we have introduced is a direct sum of a massless, spin-1 unitary representation (corresponding to 
the indices 1 and 2) and a massless spin-0 unitary representation (corresponding to the indices 3 and 0) of the Poincare 
group. In such a structure the 'longitudinal' and 'timelike' components do not transform as parts of a four-vector, but 
as two scalar fields. However, there is a lot of freedom here. There exists, for example, an interesting representation 
which employs the link between Bogoliubov and SO(l, 3) transformations, and where all the four components behave 
as those of a four- vector. A particular example of this link was employed in Sec. V. It is interesting to compare the 
two representations in the context of field quantization, but this will be done in a separate paper |3J|. 

XI. LORENZ CONDITION AND EUCLIDEAN GROUP 

The field tensor F^^x) = d a A b {x) — dbA a {x) consists of two parts corresponding to spin-1 and spin-0 fields (formulas 
(JEIi and JE3J, respectively) 

Fab(x) = -JdT(k)Tr A (k)ir B (k)e A >B> (a(fe, -Je"**'* + a(k, +) V*'*) + h.c. (83) 
- V2/dr(fc)*Af afc (fc)(n i (fe)cosfc-a; + n 2 (fc)sinfc-a;) (84) 

*M ab {k) = -k a uj b (k) + k b uj a (k) (85) 
n x (fe) = i(a(fc,3) + a(fc,3)t + a(fc,0) + a(fc,0)t) 
n 2 (fc) = ^(a(fe,3) -a(fc,3) f + a(fe,0) t - a(fe,0)) 

The tensor * M ab (k) is the dual of 

M ab (k) = iir (A LL> B )e A >B> ~ m{A'^B^AB (86) 
k b *M ab (k) = -k a (87) 

One immediately recognizes in (18611 and (|87|) spinor formulas for a massless angular momentum tensor and the Pauli- 
Lubanski vector of helicity —1 (cf. [13, Eq. (6.3.2)). 

The gauge transformation Q influences the part (|84|l in F_ ab (x) according to 

*M ab (k) h+ * M ab (k) - k a (k)q b (k) + k b {k)q a {k) 
q a (k) = <p(k)m a (k) + 4>(k)m a {k) 

that is, in a way typical of angular momentum. The 4-vector q a {k) can be used to reexpress the gauge transformed 
spin-frame as a twistor 32] 

n A (k) = n A {k) (90) 

u A {k) = CJ A {k)+q AA ,(k)TT A '(k). (91) 
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Mutual relations within an equivalence class are thus determined by the twistor equation. Change of origin in the 
space of coordinates q a can be compensated by the Bogoliubov transformation B_. 
Together with J_(k) occuring in (|75|l we obtain the algebra e(2) 

= iS r (k,k')Il 2 (k) (92) 
[J(k),U 2 (k')} = ^(fe.feOn^fe) (93) 

EI 1 (fc),n a (fe')] = o (94) 

It is interesting that the removal of the scalar fields by the constraint 

(^|n 1 (fc)|*) = (*'|I! 3 (fc)|*)=0 (95) 

is analogous to the condition leading to the classical Maxwell field if one starts from induced representations [3l| . 
Indeed, all massless discrete-spin representations are found if one requires that the two translation generators of e(2) 
annihilate vectors from the representation space. It might be therefore tempting to impose the stronger constraint 

n 1 (fc)|*)=n 3 (fc)|*)=0 (96) 

also here. To see why this condition would be too strong we write the potential in terms of e(2): 



A a (x) = -V2jdT(k)uj a (k)(ll 2 (k)coskx- n^fe) sinfcx) - \ / 2jdT(k)k a (Q 1 (k) coskx + QJk) sm.kx) +...(97) 
where 

Qi( fc ) = -5?(«(fe J 3)-a(fc,3) t -a(fc,0) t +a(fc,0)) 
Q 2 (fc) = -i(a(fe,3) + a(fe,3)t-a(fe,0)t-a(fe,0)) 

and the dots stand for the part involving only the spin-1 fields. The part involving Q(k) is a gauge term and this is 
why we do not see it in F_ ab (x). The entire algebra reads 





= i5 r (k,k')E 2 (k) 


(98) 


[J(fc),n 2 (fe')] 


= i<5 r (fe,fc')ni(fc) 


(99) 


[Q 1 (fc),-J(fc')] 


= ^ ( 5^(fe,fc')Q 2 (fe) 


(100) 


[-j(fc),0 2 (fe')] 


= i* r (fc,fc / )Q 1 (fc) 


(101) 


[Q i (fc),n 1 (fc')] 


= iS T (k,k')I(k) 


(102) 


[Q 2 (k),Ua(k')} 


= i5 r (k,k')I(k) 


(103) 


[n 1 (fe),n 2 (fc')] 


= 


(104) 


[Q^klQ^k')} 


= 


(105) 


[Q 1 (fc),2 2 (fc')] 


= 


(106) 


[Q 2 (k)Mk')} 


= 


(107) 



This is the Lie algebra of the 2-dimensional Euclidean group in a phase space. The "position operators" Q(k) shift 
the "momenta" II (fc) and the constraint (|96|l must be inconsistent with dynamics. There is no problem with (|95|l . 
The 4-divergence of the potential 

d a A a (x) = \/2/dr(fc)(n 2 (fc)sinfca; + n l (fc)cosfca;) 
shows that the weak Lorenz condition 

(^'\d a A a (x)\'f) = (108) 



is equivalent to 19511 . 
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XII. LORENZ CONDITION AND POINCARE COVARIANCE OF STATES 

The 1-oscillator Hilbert space H consists of vectors 

I*) = EZ,n u n 2 ,n 3 =oI dT ( k )^( k ^ n o^ n ^ n 2,ri3)\k,n a ,n 1 ,n2,n 3 } (109) 

satisfying 

En , ni ,„ 2 ,„ 3 =o/^( fe )l*( fe ; «o. ni.na, n 3 )\ 2 < oo (110) 
Subspaces consisting of vectors of the form 

|*no,rn,n2,n 3 ) = I^(k)^(k, n , n x , n 2 , n 3 )|fc, n , ni, n 2 , n 3 ) 

are invariant subspaces of the representation constructed in Section Ixl 
In particular, all the vectors of the form 

I*) = Er i ,n 2 =o/ tff ( fe )*( fe ' n o>«i>«2,n3)|fe,no,ni,n2,n 3 ) (111) 

belong to the Poincare-invariant subspace satisfying the weak Lorenz condition (|95|l for N = 1. Moreover, if we 
additionally require no = n 3 then 

(*'|Pf|f ) = (112) 

An extension to arbitrary N is immediate. One concludes that the Lorenz condition Q95JI can be imposed in a Poincare 
invariant way. 

XIII. VACUUM, MULTIPHOTON, AND COHERENT STATES 

The subspace corresponding to no = n\ = n 2 = U 3 = defines the vacuum for N = 1. Any vector of the form 

|0) = /dT(fe)O(fe)|fe,0) (113) 
plays a role of a 1-oscillator vacuum. For arbitrary N the vacuum state is taken in the form 

\0) = \0)®---®\0) (114) 

S v ' 

AT 

All vacuum states are annihilated by all annihilation operators. Vacuum states are translation invariant and SL(2,C) 
covariant: 

UaJO) = Jrfr(fc)O(A- 1 fc)|fc,0) (115) 
HaJQ) = U A ,y\0)®---®U A ,y\0) (116) 

Of particular importance is the scalar held representing vacuum probability density Z(k) = |0(fc)| 2 . Square integra- 
bility implies that Z(k) decays at infinity; later on, we will also require Z(k) going to zero at k — in order to avoid 
infrared divergences. The latter would spoil the explicit constructive nature of the present approach. The number 
Z = max/j{Z(fc)} is Poincare invariant and can be interpreted as a renormalization constant. 

Multiphoton states are obtained in the usual way by acting on the vacuum \0) with creation operators. Coherent 
states associated with amplitudes a(k, ±) occurring in (0) are defined in terms of the displacement operator [|| 

V(a) = cxp (a(a) f - a(a)) (117) 

«(") = E*=± J^T (k)a{k, s)a{k, s) (118) 

A coherent state is constructed from vacuum by \0(a)) — r D_(a)\Q). Coherent-state averages are related to classical 
fields by 

(Q(a)\A a (x)\0(a)) = iJdT(k)Z(k)(m a (k)a(k, +) + m a (k)a(k, -))e' lk - x + c.c. (119) 
Let us note that the averages involve the amplitudes Z(k)a(k, ±) and not just a(k, ±). 
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XIV. FIELDS PRODUCED BY A CLASSICAL CURRENT 



The interaction Hamiltonian in the interaction picture is assumed in the usual form 

H(t) = Jd 3 xJ a {t,x)A a {t,x) (120) 

where J a (t, x) is a classical conserved current. The interaction picture evolution operator satisfies 

i£ t U(t,t ) = H(t)U(t,t ), U(t ,t )=I (121) 

Recalling that A a (i,:r) depends on time via the free Hamiltonian Hq = P (|72() we can split the Heisenberg-picture 
time evolution into parts involving separately the interaction picture U(t,to) and the free evolution, i.e. 

A?{x) = U(t, t )^A a (x)U(t, t ) (122) 

To obtain the latter we made the usual assumption that there exists a time to at which the field is free. This restriction 
is eased later on by moving to to ±oo. The following two splittings of (|120|) are important 



H(t) = H x (t) + ffi(t)t = H 2 (t) + H 2 (tf 

H x {t) = iJd 3 xJ a (x)JdT(k)(x a (k)a{k,l) + y a (k)a(k,2) + z a {k)a{k,3) + t a (k)a(^ (123) 

H 2 (t) = ifd 3 xJ a (x)jdr(k)(x a (k)a(k,l) + y a (k)a(k,2) + z a (k)a(k,3)-t a (k)a(k,0)e 2th (124) 

since the commutators 

[ff<(ti),ffi(*2)] = (125) 

[fl-i(ti)t,fli(t 2 )t] = (126) 

[Hifa^Hjfo^Hjfoy]] = (127) 

[H^WH^.H^] = (128) 

hold for all i,j = 1, 2 and arbitrary times. The commutators 

[i?i(x ), J ffi(yo) t ] - ifd 3 xd 3 yJa(x)D ( - +) (x-y)J a (y) (129) 



[H 2 (x ),H 2 (y ^] = ifd 3 xd 3 yJ a (x)D(+\x~y)J a (y) 

+2 J d 3 xd 3 y J a {x)J b (y) J dT(k)I_(k)t a (k)t b (k) cos k ■ (x - y), (130) 

are in the center of CCR. Employing continuous Baker-Hausdorff formulas |3^, 03 

Texp (f to dr{A(r) + B(r))) = exp (^ Q dr(A(r) + B(r))) exp (§ /'/n J^dr, (6(n r 2 ) - 9{r 2 n )) [A( n ), B(r 2 )}) 

= exp U\drA{r)) exp U\ o drB{r)) exp {^dr^dr^ - T 2 )[B(n), A(r 2 )}) (131) 

where A(r), B(t) satisfy relations analogous to (|125|I - H128|I . we find that 

U(t,t Q ) = exp(- l /^ 4 xJ Q (xUJx))exp(-i/^^Vd 4 x 2 J a (xi)^ adv (x 1 -x 2 )J a ( a ; 2 )), (132) 

where i2adv( x ) = ^S(—xo)D.{x). Formula 1)131(1 will be later used to compute the photon statistics. 
Employing 1)1321) we find the explicit form of the Heisenberg-picture evolution 

A?{x) = Mx) + Jl^yRix - y)J a (y) (133) 

Our field Af (x) is free at t = to- In the next Section we show that the weak Lorenz condition holds for (|133fl only in 
the limit to = ±oo. 
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XV. LORENZ CONDITION AND RETARDED/ ADVANCED SOLUTIONS 

The four divergence of (|133ll takes the form 

d a Aa( x ) = free part + Jd s x'D(t - t ,x - x')J (t ,x') 

Taking an arbitrary coherent-state average 

(0(a) | 0M? (x) 1 0(a)) = Jd 3 x' (0\D(t - t ,x - x') \ O) J (t ,x') 

and requiring the Lorenz gauge for all conserved currents, we obtain a condition on the vacuum-state probability 
density Z(k) = |0(fe)| 2 

(Q\D(t - t , x)\0) = i f dT(k)Z(k)(e-' k - x e lkoto - e - ik ° to e ik - x ) = 
This cannot hold in general if to is finite. However, for to — > ±00 the condition becomes equivalent to 

lim J d 3 kf(k)e i ^ to = (134) 

where f(k) = Z(k)e lk ' x /\k\. The latter condition requires only that Z(k)/\k\ satisfies assumptions of the Riemann- 
Lebesgue lemma. 

We thus restrict the analysis to the two cases of either retarded or advanced solutions. The formulas are 

(x) = A a (x) + Jd^yD^x - y)J a (y) (135) 

At v (x) - A a (x) + /d 4 y^ adv (x - y)J a (y) (136) 

D Iet (x) = 8(x )D(x) (137) 

^adv(^) = ~0(-x o )D(x) (138) 

D(x) = D rct (x)-D &dv (x) (139) 



Since DD_(x) = we find 



OD Iet (x -y)= DD adv (x - y) ^ S(x - y) (140) 



One has to bear in mind that S(x — y) is defined by H14()|) and that the resulting operator is not equivalent to the 
Dirac delta. 

The advanced and retarded potentials satisfy 

□Af /adv (x) = Jd±y5(x - y)J a (y) d = f J» (141) 
The weak Lorenz condition implies that the average current 

(J a (x)) = (0(a)|J a (z)|0(a)) = (0\J^(x)\Q) (142) 
is the conserved physical current that produces the classical electromagnetic field 

(A ret/adv (x)) = (0(a)\AT^ v (x)\0(a)) . (143) 

The modification of the current depends only on the choice of the vacuum state because the displacement operator 
commutes with £(fc). 

To close this Section let us mention that an operator analogue of the Feynman propagator 

D F (x) = 9(x )D (+) {x) - 6(-x )D(-\x) (144) 
= R adv (x) +#(!) =R Iet (x) -rt-\x) (145) 

would occur in perturbative formulas in exactly the same places as in the standard formalism. The reason is that the 
algebraic structure of Feynman diagrams is unchanged by the change of representation of CCR. Since □_D^ ± - ) (a;) = 0, 
the Feynman potential 

Al(x) = A a (x)+Jd 4 yD F (x-y)J a (y) (146) 
satisfies the same equation as the retarded and advanved fields, but is non-Hermitian for real currents. 
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XVI. DYNAMICS OF RETARDED AND ADVANCED SOLUTIONS BETWEEN TWO FINITE TIMES 

We have solved the Heisenberg equations with free-field "initial" conditions at to = ±00 and arrived at retarded 
and advanced solutions of Maxwell's equations. We have not yet shown what kind of dynamics will map retarded 
or advanced solutions at a finite time t\ into retarded or advanced solutions at another finite time t. This would be 
the true solution of the Heisenberg-picture evolution since at a finite initial time the field cannot be free, unless the 
charge of the current is zero. 

One can immediately write down appropriate formulas on the basis of the retarded and advanced solutions 

AT^ dv (t,x) = lira U^toYUo^hyuit^t^AT/^it^x^it^toYUo^t^U^to) (147) 

to — > — ocf +00 

= W_ /+ (t,t 1 YA! a ct / adv (t 1 ,x)W^ /+ (t,t 1 ) (148) 

Let us recall that Uo(t,ti) = cxp ( — iHo(t — £1)), where Ho = P ^ s the free Hamiltonian defined by the reducible 
representation of CCR, and U(t, to) is the interaction-picture evolution operator. Some care is needed in the definitions 
of W±(t,ti) if the limits lim to ^± oa U(ti,toyUo(t,ti)U(t,to) involve divergent phase factors. This is the standard 
problem and has nothing to do with the divergences of quantum field theory. Keeping this subtlety in mind we arrive 
at 



W±(t,ti) = exp(i/* ± 1 co d 4 a;(j a (xo,a;)- J a (.T +t-t 1; a;))A a (a; ,a;))exp(- 2 Ho(t-<i)) (149) 

It is clear that for a static charge density the evolution is free. The ranges of integration are finite also in case the 
currents are static for t < t_ and t > t + with some t± . It should be stressed that this type of "switching on and off" 
of the current is perfectly consistent with charge conservation. 
The corresponding Hamiltonian H± (t) satisfying 

id t W±(t,ti) = W±(t,h)H±(t) (150) 

reads 

H±(t) = H + f t ±oo <I l x£ a (x) 1 g s J a (x). (151) 

In the next section we show the explicit form of H± (t) for a pointlike charge. As we shall see the Hamiltonian has a 
clear physical interpretation. 



XVII. EXPLICIT FORM OF THE NEW HAMILTONIAN FOR A POINTLIKE CHARGE 

A pointlike charge q localized on an infinitely long world- line z a (t) = (t, z(t)) leads to the conserved current [35ll3r| 

J a (t,x) = q(l,v(t))S^(x-z[t)) (152) 

where v(t) — dz(t)/dt. Let us assume that the world-line represents a charge which is at rest for times t < t_ and 
t > t+. The assumption implies also that v(t±) = if we assume that t 1— > z(t) is twice differentiable. Under these 
assumptions we find for t- < t < t + that 

H± = Ho- qA(z(t)) ■ v(t) - qf^dz ■ E. (153) 
H- = Ho for t < t- (154) 
H+ = H for t > t+ (155) 

The line integral in the third term of (|153l) is along the part of the charge world-line where the charge velocity is 
nonzero. The electric field operator takes the usual form E_ = — do A — VAq. 

Let us note that the terms explicitly involving A_o have cancelled out. It is clear from the construction that the 
electric field occuring in H± is free. Therefore, the Hamiltonian does not contain self-energy terms but, instead, takes 
into account the work performed by the particle against the electric field. 
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XVIII. PHOTON STATISTICS 



The operator U(t, ±00) (as well as Uo(t,±oo)) in general does not exist due to the problem with divergent phase 
factor. Fortunately we do not really need U(t,to) itself, but only its action on operators X 

U(t,t )(X) = U(t,t )iXU(t,t ) (156) 

Similarly, in order to compute the S-matrix we concentrate on the limiting operator map S — U(+oo, —00). 
Eqs. {T55Jl. JHll, CCS>t imply that at one hand 

if (x) = lim lim U(r,T )^ dv (x)U(r,r ) = S(A^ v (x)) (157) 

Tq — >— OO T — »- + 00 

and on the other 

A a {x)+Jd i yD{x-y)J a {y) = S^x)). (158) 

Finally, as shown in 0] , the S'-matrix S gives the action of the displacement operator on the field operators. 

More interesting is the question of photon statistics in fields produced by classical currents, especially if the currents 
are pointlike or stationary. Accelerated pointlike charges lead, in the standard formalism, to infrared catastrophe. 
In a naive approach, all transition probabilities are zero, which contradicts unitarity of the ^-matrix. In manifestly 
covariant approaches, such as the Gupta-Bleuler formulation, static charge distributions lead to infinite vacuum-to- 
vacuum probability, which again makes no sense. A mathematically correct treatment is possible [2^, but requires 
a change to an abstract representation. This goes against the philosophy of explicit construction, adhered to in the 
present paper. 

In our approach infrared divergence can be avoided by requiring that the vacuum probability density Z(k) tends 
to zero at k = 0. Let us take the operator H 2 (t) given by (|124J) and split it into the parts H 2 30 \t) and H^ 2 \t) 
involving, respectively, the fields of spin-0 (i.e. a(fc,3) and a(fc,0)) and spin-1 (i.e. a(k, 1) and a(fc,2)). Let us recall 
that the ./V-oscillator Hilbert space is spanned by iV-fold tensor products of vectors of the form |fe, ni, 712,713, no). We 
shall refer to such 1-oscillator states as containing rii + n 2 transverse excitations and n 3 + n longitudinal ones. Any 
state belonging to the iV-oscillator Hilbert space 7t and containing n transverse excitations, where n is is the sum of 
the transverse excitations of all the ./V oscillators, is regarded as a state involving n transverse photons. Similarly we 
define a general state involving n' longitudinal photons. In particular, the vector 

H { 2 12 \t^ ...H { 2 12 \t n )^\0) (159) 

belongs to the subspace of n-transverse-photon states. The state 

H^ihY ...H { 2 30) (t n ,y\O) (160) 

involves n' longitudinal photons. 

Denote by P n ri the projector on the subspace of H that contains states with n transverse and n' longitudinal 
photons. The probability of finding n transverse and n' longitudinal photons in the state produced from vacuum by 
a classical current is thus 

ftm'(Mo) = {O\U{t,0Vnn'U(t,t O )\O). (161) 

Employing (|131J) we find for p rm > — p n n'(oo, —00) 

I / 1 H n d n ' 

nln'r— 1 12 30 '— ' nln'l dfj, n dv n v ; 

C(p,v) = (O\e» Fl2 e vF30 \O) 

Assuming that detectors react only to spin-1 photons we obtain photon statistics 

00 

Pn = V Pnw = ~. (Q\F? 2 e" Fl2 \Q) (162) 
z — ' n\ 

n'=0 

Alternatively, in order to describe quantum optics of spin-1 observables directly, without any reference to the spin-0 
fields, we can consider states defined via reduced density matrices with the spin-0 parts traced out. Probability p n is 
an example of an average computed in terms of such a reduced density matrix. 
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A. Fourier description 

To find an explicit formula we first consider a current whose 4-dimensional Fourier transform is a well behaved 
function J a (k) = J d 4 xJ a (x)e ik - x . Then 

F 12 = J dT(k)I_(k)(\J a (k)m a {k)\ 2 + \J a {k)m a {k)\ 2 ) 

F 30 = J dT{k)I_(k)\J a {k)u a {k)\ 2 

where J a (k) is a restriction of J a (k) to the light-cone. Due to the continuity equation J a (k)k a (k) = the spin-1 
expression F\2 is independent of the choice of LOA{k). 

Employing the relation between 7(fe) and I{k) = \k){k\ ® 1 we can write the generating function as 



C(jm,u) = \^dr(k)Z(k) C x V ^(\J a (k)m a (k)\ 2 + \J a (k)^ (163) 

Of particular interest is the thermodynamic limit N — > oo for the spin-1 part. The corresponding generating function 
becomes 

lim^C(/i,0) =c^^JdT(k)Z(k)(\J a (k)m a (k)\ 2 + \J a (k)fh a (k)\ 2 )^) (164) 

and 

p n = ±^dT(k)Z{k)(\J a {k)m a {k)\ 2 + |J Q (fc)m a (fc)| 2 ))"cxp ( - JdT(k)Z(k)(\ J a {k)m a (k)\ 2 + \J° \k)m a {k)\ 2 )) 

This is basically the well known Poisson distribution, with one modification: The standard infrared-divergent result 
is found if one puts Z(k) = 1. However, we know that JdT(k)Z(k) = 1 and thus Z(k) ^ 1. We have remarked 
earlier that the maximum value of Z{k) is a positive Poincare invariant, denoted by Z . Introducing the new function 
x(fc) = Z{k)/Z, and absorbing Z 1 / 2 into a renormalized current J r " n (fe) = Z 1 / 2 J a {k) we find 



Pn - ±(jdT(k) X (k)(\J? cn m a (k)\ 2 + |J r a cn m Q (fc)| 2 ))"cxp ( - JdT(k) X (k)(\J^ n m a (k)\ 2 + \J? cn rh a (k)\ 2 )) 

Now this is indeed the standard regularized expression. The latter provides us with a new information about the 
vacuum wave function 0(k): It has to vanish at the origin k = if one wants the cut-off function x(&) to regularize 
the infrared divergence. The origin belongs to the boundary of the light cone. Vanishing at the origin is a Poincare 
invariant boundary condition. 

We regard this result as very important, as it handles in a natural manner two elements that are imposed in an 
ad hoc manner in the standard formalism. First of all, we do not need to justify the infrared cut-off by hand-waving 
arguments on unobservability of "soft photons". Our formalism introduces the cutting-off function automatically. 
Secondly, we know what is the origin of a renormalization constant: This is simply the Poincare invariant associated 
with the vacuum wave function. 

B. Pointlike static charge 

In this case it makes no sense to switch to the Fourier domain, since the position space calculation is more straight- 
forward. Assume the current is J a {x) — (qS^(x),0). The generating function becomes 

cm = (jdr(fc)z(A^ 2 ^^ (165) 

In the thermodynamic limit 

lim C(ji,v) = exp(/ ( ir(fc)Z(fe)(2g 2 /|fc| 2 )( Ma; o(fc) 2 + Myo(fc) 2 + ^o(fc) 2 + ^o(fc) 2 )) (166) 
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Identical results are obtained if instead of U(oo, — oo) one works with U(t, ±00) for a finite t. 

Let us remark that an analogous calculation performed in a reducible version of Gupta-Bleuler formalism [17| leads 
to vacuum-to- vacuum "probabilities" that are greater than 1. The reason is that for currents whose only nonzero 
component is Jo(x) the Fourier-space version of continuity equation does not read J a (k)k a — 0, but Jo(fco)^(fco) = 0, 
and one cannot claim that J a {k) is spacelike. In our formalism the timelike component of the current comes with the 
correct sign. 

C. Renyi statistics for finite N 

Generating functions can be written in a unified way for any N in terms of Kolmogorov-Nagumo averages of the 
form used in Renyi statistics. Let us recall that Renyi's alpha entropies were obtained in |37| as Kolmogorov-Nagumo 
averages 

of the random variable Ij = ln(l/p J ), and (f>(x) — e^ a ^ x . For a — 1 one obtains the standard Boltzmann-Shannon 
entropy. In |38j ] it was shown that thermodynamics that employs Renyi type averaging can be used to derive certain 
equilibrium distributions occuring in linguistics and protein folding. Various arguments based on thermodynamics 
suggest that a / 1 statistics may be typical of finite systems. Photon statistics for finite-iV representations of CCR 
supports this intuition. 

Indeed, in the thermodynamic limit we found a generating function of the form C(/i, 0) = e" W' with 

= fiJdT(k)Z(k)(\J a (k)m a (k)\ 2 + \J a (k)m a {k)\ 2 ) 

being a linear (a = 1) average of fi(j J a (k)m a (k)\ 2 + |J a (fc)™ a (fc)| 2 ) , with probability density Z(k). For finite N we 
find C(n,0) = e where <j)(x) = e^ a ^ x , a = 1 — 1/N, and we average the same random variable with the same 

probability distribution. Obviously, the limits N — > 00 and a — ► 1 are equivalent. It follows that the field theories 
based on N < 00 or N = 00 representations are related to one another in a way that is analogous to the relation 
between systems described by < a < 1 and a = 1 entropies. These, on the other hand, are known to apply to 
fractal and non-fractal geometries, respectively. A natural intuition thus relates the N < 00 case to some "space-time 
foam" , and N = 00 to continuum space-time. 

XIX. CLASSICAL FIELDS PRODUCED BY CLASSICAL SOURCES: A QUANTUM WAY 

Our previous analysis shows that, having a classical current J a (x), we obtain a result that agrees with standard 
calculations if one (a) absorbs Z 1 / 2 in the current by means of J™ n (x) = Z 1 / 2 J a {x) (bare charge renormalization 
q 1 ► q Icn = Z 1 / 2 q), and (b) compares the result with a regularized formula, which is anyway the one we have to 
compare with experiment. Z is not a constant but rather an invariant of the Poincare group that characterizes 
a given vacuum. We have also obtained a cut-off function x{k) = Z(k)/Z. At this stage we do not have much 
information as to the exact form of and can only say that it vanishes for large fc and k — 0, and that x(fe)/|fe| 
fulfills the assumptions of the Ricmann-Lebesgue lemma. 

Now let us take an arbitrary classical amplitude a(k, ±) corresponding to left- and right-handed Fourier modes of 
a classical electromagnetic field. We define the quantum optics regime by the support of those classical amplitudes 
that satisfy 

a(k,±) = a(k,±) X (k). (168) 

The latter formula is meaningful provided = 1 if fc belongs to quantum optics regime. Classical fields belonging 
to quantum optics regime do not contain wavelenghts that are either too large or too small. Let |0(a)) be a coherent 
state with a in quantum optics regime. Let us take the coherent state average of the retarded solution of Heisenberg 
equation of motion i|135[l and express it in terms of the renormalized current J™ n (x). Taking into account that 



A a {x) = (0{a)\A a (x)\0(a))^ZiJdT(k)(m a (k)a{k,+) + m a {k)a(k,-))e- lk - x +cx. 1 (169) 
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by i|119|) and the assumption that i|168[) is fulfilled, the formula 



(Q(a)\D(x - y)\0{a)) = Zifdr(k)x(k)(e~ lk - (x ~ y) - e^*"^), (170) 
and dividing the entire solution by Z 1 ' 2 , we find that the classical field 

(A?(x)) = Z- 1 / 2 A a (x) + Jd 4 yD Ict (x~y)JS h y s (y) (171) 

d a (Al ct (x)) = d a A a (x) = (172) 

d a J% hys (x) = (173) 



exhibits the textbook relation between the in-field, renormalized current, and renormalization constant. 

Here D rct (x — y) is the ordinary retarded Green function and Ja hys (y) is the efective current obtained after charge 
renormalization and inclusion of x(&) i n its Fourier transform (the convolution of D r „ t (x — y) and J a {y) in i|135|) allows 
to shift the regularization from the Green function to the current, and vice versa). All these objects have occured in 
our calculation automatically. 

Finally, let us have a closer look at the effective current. For simplicity take a static pointlike charge. Employing 
the relations 

Z 1 ' 2 Jt ys (x ,x) = J d 4 y(0\5(* - y)\Q)Jo(y) = q(Q\d D(0, x)\0 = Zq\ J j0zx(k){e^ s + e"**) 
we find that the effective total charge 

Q = | d^xJ^ixo-x) = Z^ 2 q X (0) = <f cn X (0) 

vanishes since we require x(0) = 0- 

In order to check the physical meaning of this condition let us consider the simple case of spherically symmetric 

XkiteK*) - | o otherwise { ' > 

Then Q = q Ten for k\ = and Q — for hi > 0. One can check that the effective charge density ^^(xq, x) = 
Pk 2 ( x ) — Pki( x )j where J d 3 xpk 2 (x) = J d 3 xpk 1 {x) = q lcn for all k2,ki > 0, but simultaneously the pointlike limit 
limfe^o Pk t i x ) = holds. It turns out that the charge density consists of a difference of two densities: One, which is 
the sharper and more localized the greater &2, and the other which is the flatter and less localized the smaller k\. One 
of them corresponds to localization of the charge g rcn in a sinc-like very sharp-peaked density, and the other describes 
the charge —q rcn distributed almost uniformly in a volume which becomes infinite if fci = 0. Plots of the densities 
for ki > and k\ = become indistinguishable even for relatively large fci and small k2, so we leave this exercise to 
the readers. In practice, one cannot locally distinguish between k\ = and k\ w 0, but globally the two cases are 
inequivalent. 

XX. SUMMARY AND CONCLUSIONS 

We have discussed a new quantization scheme based on reducible representations of CCR. The principal goal of this 
research program is to arrive at a mathematically consistent formalism for quantum fields, that should be described 
by a theory and not a set of working rules, as Dirac summarized the current status of field quantization |3j]. In our 
approach fields are represented by operators and not operator- valued distributions. The field is a finite quantum 
system, and the measure of its size is the parameter N. For this reason two sources of infinities are absent in the 
formalism from the very outset: We can multiply field operators at the same points in space-time and all the tensor 
products one encounters are finite. The latter condition means that we deal only with factors of type I, in von 
Neumann's terminology. 

We have carefully analyzed Poincare covariance of the theory. There were two aspects we had to understand to 
make sure that the new quantization is not inconsistent with special relativity. First of all, we constructed a unitary 
representation of the Poincare group whose carrier space is an ordinary Hilbert space involving no indefinite metric. 
Four-potential operator is self-adjoint, the dynamics is unitary, but commutation relations for fields are nevertheless 
manifestly Poincare covariant. The formalism is a promising alterantive to the Gupta-Bleuler quantization, where the 
price payed for manifest covariance is either in non-positivity of the Hilbert-space metric, or in non-Hermicity of the 
potential. 
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Secondly, we had to understand in what way a Lorentz transformation influences the gauge freedom. The latter 
has led to the observation that a change of gauge due to Lorentz transformations can be always compensated by a 
Bogoliubov unitary transformation of the vacuum. An inclusion of the Bogoliubov transformation turns the 4-potential 
into a 4-vector field. 

One element that remained arbitrary is what kind of a tetrad one has to associate with a 4-potential. Once one 
makes a choice then the remaining freedom can be controlled by Bogoliubov transformations, whose explicit form has 
been given. In our formalism the dynamics is unitary, in the ordinary meaning of this word. We have no problems 
with negative or greater than 1 "probabilities" that occur in the Gupta-Bleuler formalism. The correct probability 
interpretation is guaranteed by the Schwartz inequality. Further, if one looks at the radiation fields then our formalism 
produces the standard regularized formula, which does not depend on a choice of gauge. 

Finally, its seems that we have produced the first example of Heisenberg dynamics where the retarded or advanced 
fields unitarily evolve from, say, t = to another finite t. Our construction allowed to systematically derive the 
form of Hamiltonian that is responsible for such an evolution, and the result turned out to differ from the usual 
minimal-coupling expression: There is no scalar-potential part and a new term occurs. The term describes the work 
performed by the charge moving in electric field. This result may have implications for quantum optics where the 
usual treatments of spontaneous emission or resonance fluorescence are based on initial conditions at t = and not 
t = -oo (cf. [UEll], the exception is H2). 

Summing up, we think we have proposed at least a nontrivial answer to the problem posed by Dirac in his last two 
papers. It looks like the formalism, supplemented by its fermionic analogue introduced in is ready for calculations 
in full quantum electrodynamics. Some preliminary results on loop diagrams have been already obtained and will be 
reported in a future paper. 
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